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Abstract 


We propose a parametrization for interpreting some of the presently available 
data of the —)• K^pp decay, in particular those by LHCb and Belle 

collaborations. The model is inspired by the well-known current and transition 
contributions, usually assumed in this kind of decay. However, in the light of 
considerations about the dominant diagrams and about final state interactions, we 
modify some parameters of the model, determining them by means of a best fit to 
data. We show the results, which we discuss in some detail. Moreover we give some 
predictions on other observables relative to the decays. 
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1 Introduction 


The physics of the 5-meson has opened a new door in the sector of hadronic weak 
interactions. In particular, numerous B decays, analyzed during more than a decade, 
have confirmed the CKM mechanism for CP violation[T]; yet, it is generally believed that 
some particular decays of this kind could reveal new physics beyond the Standard Model 
(SM). Therefore several experiments in this sense have been suggested or even realized. 
In particular, data of —)■ K^pp decay have been published recently by the LHCb 

collaboration [2-4]. 

Three-body baryon-antibaryon B decays were detected for the hrst time by ARGUS 
in 1987 15]. Although immediately ruled out by CLEO|6], that claim gave rise to a lively 
theoretical interest in the subject for several years: see ref. [7] for a complete and thorough 
review. More recently, this interest was revived[8-12] and new experimental results were 
found by the CLEO[l3], Belle[Tll[T5] and BaBar[T6l HTj collaborations. 

A three-body decay presents some advantages over a two-body one. For example, if 
one chooses a three-body charged B decay such that all final particles are charged, its 
detection is favored. Moreover, an amplitude analysis in different intervals of the Dalitz 
plot allows direct measurement of the relative strong phases [T^. 

Baryonic three-body decays of charged B are especially useful for studying the strong 
dynamics and, in principle, for analyzing the observables sensitive to new physics[8]. 
Indeed, it is customary to assume that hnal state interactions [I9] (ESI) are much smaller |2j 
than those of three-meson decays, characterized by large re-scattering effects [T9l 1201 EH]- 
This fact would allow the study of ESI by comparison between the two different kinds 
of decays [2] and therefore could help in determining the weak phase, which generally 
presents serious difficulties [21]. However, as we shall see in the present paper, ESI might 
play a surprisingly important role also in baryonic decays. 

Besides, such decays are somewhat advantageous with respect to the baryonic two- 
body ones (e. g., 5° —)■ pp), in that they give rise to greater decay rates than two-body 
baryonic decays. This is connected with the threshold baryon-antibaryon peaks, predicted 
by Hou and Soni[8| as a consequence of the increasing difficulty of a quark to hadronize to a 
baryon at higher energies. In particular, if the final state includes the proton-antiproton 
system, the effect is attributed to re-scattering[22|. Moreover, the decays considered 
exhibit intriguing forward-backward (FB) and Dalitz plot asymmetries [T71IT^ 1^ . 

This kind of decays is generally described by assuming two factorizable amplitudes of 
the type PE] 

r= {M\{qsq2m{BM{qib)\B-) and I = {B,B2\{qiq2m {M\{q-,b)\B-), 

called respectively transition and current terms. Here M, Bi and B 2 are the final meson, 
baryon and antibaryon respectively, while b and q {q) are the destruction (creation) 
operators of the active (light) quarks involved in the decay. 

Turning to the decay B^ —)■ K^pp (to be named K decay in the following), it presents 
two main differences with respect to the B^ —)■ ir^pp (tt decay from now on). Firstly, it 
has a wider pp peak, which may be explained as an effect of the interference between the 


1 


transition and the current contributions |9] : indeed, while suppressed in the vr decay, the 
current term contributes considerably to the K decay and decreases more slowly than the 
transition term at increasing pp effective mass. Secondly, the sign of the FB asymmetry 
{Aps) in the vr decay is opposite to that in the K decay: while the behavior of the vr 
asymmetry can be explained in the framework of quark dynamics[7j, the other one is more 
difficult to interpret. 

An important feature of the K decay is that it derives contributions from two different 
amplitudes, typically tree and penguin, endowed with different weak and strong phases. 
Therefore, we expect that data exhibit a sizable direct CP asymmetry (Acp), to be 
compared with the SM predictions. Indeed, this observable has been measured recently to 
high precision by the LHCb collaboration [2]. Unfortunately, as recalled above, comparison 
with theory is particularly difficult in weak decays, owing to the problem of disentangling 
real CP-violation effects from strong FSI. Therefore, in the interest of understanding and 
parametrizing as precisely as possible the effects of non-perturbative QCD in hadronic 
weak decays, we perform a best fit to some available data, like the differential branching 
fraction and the asymmetries App and Acp- More precisely, on the one hand, we adopt 
for the decay amplitude the parametrization by ref. |9]; on the other hand, however, we 
modify some parameters by means of a best fit, in view of various considerations. As 
we shall see, our analysis, based on the data by Belle [U] and on the most recent data 
of LHCb [2], leads to conclusions in contrast with previous ones; in particular, we give a 
different answer to one of the main questions illustrated above, that is, the origin of the 
FB asymmetry. Lastly, we make some predictions about other observables measurable in 
principle. 

Section 2 is devoted to the definition of some observables, either presently available 
or measurable in the future, and to the theoretical function used for interpreting them. 
In sect. 3, we present our phenomenological analysis, both assuming factorization and 
relaxing this assumption; moreover we show some predictions of our model. Lastly, in 
sect. 4, we discuss the results obtained and draw some conclusions. 

2 Observables and Theoretical Function 

2.1 Observables 

We fit the data relative to the the differential branching fraction of the decay 
K^pp^ [15], i. e., 

d'B fxpp . , 

The LHCb data[2] are re-scaled in such a way that the integral over all of the ppK 
spectrum equals the total branching fraction [23]. i. e., 

Bf{B^ K^PP) = (5.9 ± 0.5) ■ 10"®. (2) 
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On the other hand, we calculate the CP asymmetry and the FB asymmetry, whose 
experimental values are, respectively [2], 



Acp = -0.022 ± 0.031 ± 0.007 

( 3 ) 

and 

Afb = 0.370 ±0.018 ±0.016. 

( 4 ) 

The latter is dehned as 

1 

1 

+ 

( 5 ) 


^ N+ + N- 


and is the number of events for which cos6p is positive (negative), 9p being the angle 
between the meson and the opposite-sign baryon in the pp rest frame. This dehnition 
holds as well for the tt decay. 


2.2 Theoretical Function 

2.2.1 Matrix Element of the Decay 

We parametrize the decay amplitude for K^pp, according to Chua et a/. IS] (see 

also refs. HH [I|), i. e., 

v< = ^(i+n (6) 

here Gf is, as usual, the Fermi constant of weak interactions and X and T are respectively 
the current and transition terms, according to the dehnition given in the introduction. 
As regards the current term, we have 

I = I 1 +I 2 , (7) 

4 

h = (A'-|i^JB-)(^ai(pp|y|0)+a,{pp|fl„|0)), (8) 

i=l 

X 2 = {\K~\Lsh\B~)aQ{pp\Rss\9i). ( 9 ) 

On the other hand, the transition term reads as 

T = {K~\Lus\Q){aT{pp\Sub\B~) -f afi{pp\Pub\B~)). (10) 

Here we have set 

L^sb = = uiPfiUL, Ll = ulPi^ul + + siP^sl, ( 11 ) 

Ll = SLlpSL, Xj = eunLlpUL + (12) 

Rfi = UR'jpUR + dRPfidu + SR'jpSR, Lgf, = sibRi ( 13 ) 

Rss SrSr^ Lug SrUr^ Sub ub, Pub u'y^b. (14) 
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Moreover the a are connected to the CKM matrix elements and to the Wilson coefficients. 
Focusing on the B~ decay, we have 

tti = VubV*^a2, 02 = -VtbV*^a3, = -VtbV^la^, (15) 

«4 = *^5 = ~VtbV^*a5, tte = (16) 

2/77'^ 

07 8 Ct^ O 4 i CIq z ' r. (1^^) 

mb{ms + TTiu) 

The numerical values of the a and of the V are listed in the Appendix. The formulae of the 
B~^ decay are obtained by taking the complex conjugates of the CKM matrix elements. 
Since Vub is the only complex CKM matrix element in the equations above, only oi, 0:7 
and 08 are involved in the CP asymmetry. 

The various (non-perturbative) matrix elements of the quark operators are connected 
to some form factors, as shown in the Appendix. As a result we get 

-^ = ^^^Ah(Pp)Ozt;(pp). (18) 

Here Pp and pp are the four-momenta of the proton and of the antiproton respectively; u 
and V are their standard Dirac spinors, normalized as uu = —vv = 2mp, where rup is the 
proton mass. Moreover, 


Oi = O 2 = ^b75, O 3 = I, 04 = 75, (19) 

with pb being the four-momentum of the decaying resonance B. Lastly, the coefficients 
[3 are 


/Si 

/S 2 

/S 3 

13a 


Here 


4 

= 2Fi ~ fK^b(y-sFv5, 

i=l 
4 

= 2Fi ^ g\pi - fKTribajFA, 
i=l 




2mr 


2=1 


fK^bOl%Fp, 


^ i=l 


m 


K 




P.,3 I 

9 a + A 
m. Am„m 


pI 


m. 


■^a)] + 


( 20 ) 

( 21 ) 

( 22 ) 


4 4 

- fxmbajFp - 2(r -h l)mpFi ^ g\pi + 2mp(rFo ^ g\pi + fxmbajFA}. (23) 

2=1 2=1 


= ±ai + a5Si^2, (i = l,2,3,4), 

l = 2pB- (l + r)(pp + pp). 


m\ — mj^ 

{Pp + Pp)^' 
3 = Pp Pp, 
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r = 


(24) 

(25) 



Pk is the four-momentum of the ii'-meson, mb(s) the masses of the quarks and rriK and 
niB the meson masses. Moreover fx is the decay constant of the i^-meson and Fq, Fi, Fa, 
Fp, Fv5, <hi, Pa and h\ {i = 1 to 4) are form factors. In particular, the - related 
to the time-like nucleon form factors, as well as the <l>j - have been set equal to zero, 
according to ref. [9]; we shall discuss this assumption in sect. 4. In this connection, it is 
worth noting that Fq^, as well as the corresponding proportionality factors 5*0,i, dehned 
in the Appendix, have the dimensions of an energy, unlike assumed in refs. UUM- 

Lastly, as we shall see in the next section, our analysis involves also the vr decay, for 
which some variations have to be performed: hrstly, one has to replace sl(r) with d^pp 
fx with /^, nis with rrid and 14s, 14s with 14d, 14d; secondly, the values of some a* and 
parameters in the form factors have to be changed, as shown in the Appendix. 


2.2.2 Calculations of Observables 


The modulus squared of the matrix element (18) reads as 


1^1" = ^ + Si243(/32ffl)|}. (26) 

z 4mp 

Here 


fill 

= ^[^PB-PpPB-Pp-ml{ml+Pp-pp)], 

(27) 


= A[2pB ■ PpPB ■ Pp + - Pp ■ Pp)], 

(28) 

^33 

= A{-ml + Pp-pp), 

(29) 

1^44 

= A{ml+pp-pp), 

(30) 

^13 

= AmpPB-{Pp-Pp), 

(31) 

f^24 

= AnipPB ■ {pp +Pp). 

(32) 


Then the differential decay width is 


dT = 


(^iA.|2n3 

'=\27i)^2E, 




i=l 


i=l 


(33) 


where Ei and Pi {i = 1 to 3) are the energies and the momenta, respectively, of the 
proton, of the antiproton and of the iL-meson in the B rest frame. By integrating over 
all variables but the momentum of the proton, we get 


dT 


ml 

4(27r)"^ms 


\M\^P, 


d^Pi 


(34) 


where ps = |p 3 |. In order to calculate the FB (or hehcity[2]) asymmetry, it is convenient 
to express energies and momenta as functions of kinematic quantities in the pp rest frame. 
To this end, we perform a Lorentz boost from that frame to the B rest frame. Moreover 
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we integrate over the azimuthal angle of the proton: since does not depend on it, 

we get 

d^= 2 % 3 ^^ ( 4 -^ 0 dTipdcosOp. (35) 

Here tt^ = TipCosOp, Tip is the modulus of the momentum of the proton and of the antiproton 
in the pp frame and 6p the helicity angle [2j, between the meson and the opposite-sign 
baryon in the same frame; moreover t = 4{mp -|- Tip) is the effective mass squared of the 
pp system, Eq = \/t + pi and 


Eo dp3 SttI 1 .^^dpl Eo 

= ++ ai- 

Taking into account energy-momentum conservation, we have 

1 


P3 = 


2m B 


[(m| - - 2t(ml - ml-) + (^''1 


(36) 


(37) 


Our analysis requires the formulae for the differential decay width and for the differential 
FB difference, as functions of the pp invariant mass mpp = y/i. To this end, we perform 
the appropriate integrations over cos6*p - respectively and ~ J_i - expressing the 
Lorentz invariant coefficients Qij in terms of the kinematic variables just dehned. The 
results are 


dr ^ dAT ^ 

dmpp ^ dmpp 


Gp mpp7ipP3 

2 A'^Ti^mB 


/(irp)|A/(irp)J. 


(38) 


Here 


3 

I{Ep) ~ ^ ^ PiGj, 

i=l 


3 

AI{7ip) = y^^piAGj, 
i=l 


(39) 


with 


pi = 8m|(|/?i|2-f (40) 

P2 = -8{m|[(mJ + 7r2)|/?i|2 + 7rJ(|/?2p]-7rJ(|/^3p + 

+ + -^pmBEo%/32^l)}, (41) 

p3 = 16mpmBP3f“^/^3fJ(/3i/?3), (42) 


Gi = 2{A+\C), G2 = \[{- + aC)ln-^^^-2C] 
3 ha I® ~ 1| 

TipB a + 1 

G3 = —r-{ 2 -aln- --), 

h I® ~ 1| 

AGi = —-B, AG2 = -2— 

h 

AG3 = -'^[{A + Ga^)ln^-^G^ + c] 

it CL 


(43) 

(44) 

(45) 

(46) 
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and 


A 

a 


B = 47r„^, C = 


Eq 

Vt' 

Eq^/I 

27lpP3 ’ 


dt 


h = 


,/ihEo^ dt^ 2t^’ 


pK 


(47) 

(48) 


Eqs. (38) allow one to calculate several observables, to be compared with data. In 
particular, we are interested in the quantities dehned in subsect. 2.1 - that is, the 
differential and total branching fractions and the overall FB and CP asymmetries - and 
in the differential CP asymmetry. 


where 7 is dehned by eq. ([^. 


•^cp — 


7-7^ 

7- + 7+’ 


(49) 


3 Phenomenological Analysis 


Here we compare the theoretical function just written with the available observables listed 
in subsect. 2.1, that is, the differential branching fraction and the overall CP and FB 
asymmetries. In particular, as regards the hrst observable, we exclude the contribution of 
the two charmonium bands, around 3.1 GeV and around 3.6 GeV. Concerning the form 


factors involved in eq. (38), hrstly we adopt the factorization assumption, like Chua et 


a/. |9]. However, as we shall see, this disagrees with data[2l[T5]. Therefore we shall propose 
a modihcation of the model, in view of some considerations about the main contributing 
graphs. 


3.1 Factorization Assumption 

As a hrst attempt, we assume factorization for the current and transition terms, according 
to refs. [3 El [m [ 12 ] and refs, therein. Indeed, this scheme generally describes the three- 
body baryonic B decays in a satisfactory way; in particular, it seemed to be supported [9| 
by the early data of the K and tt decays|25]. Therefore a comparison with the most recent 
data of those decays [21 [T5| is in order. 

In the context of factorization, the transition term has the same proportionality 
coefficient for the K decay as for the vr decay; as regards the current form factors, the 
proportionality coefficients have been calculated by Melikhov and Stech[2l] (see also ref. 
0 )- 

In the case of the vr decay, the penguin amplitude is subdominant with respect to 
the color-allowed tree diagram[T7|, even taking into account re-scattering effects to be 
discussed in subsect. 3.2. Therefore, the current term, which derives contributions almost 
exclusively from the penguin amplitude, is suppressed with respect to the transition term. 
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which consists, instead, of tree and penguin contributions. Incidentally, in that decay, also 
a subdoniinant, color-suppressed tree diagram has to be accounted for. It corresponds, in 
the case of the B~ decay, to the quark subprocess b —)■ uud, followed by recombination of 
the d-quark with the spectator u to form the 7r“, and by the fragmentation of the active 
uu pair into pp] the p is next to the vr” in rapidity space and can resonate with it as a 
A , giving a negative contribution to the FB asymmetry. 

The above evaluation about the dominating diagrams agrees qualitatively with the 
numerical results of the current form factors given in ref. [21j; indeed, the current term 
turns out to be negligibly small for the tt decayjH]. Therefore, in a simplifying assumption, 
we neglect the current term; we impose as well the functions Fvr> and Fp to vanish, 
according to ref. [9]. Then, our expression of the branching fraction depends only on 
Fa, whose proportionality coefficient is Ca', imposing the branching fraction to equal the 
experimental value |2S], 

Bf{B^ ^ Ti^pp) = (1.62 ± 0.20) • 10"®, (50) 


yields 


\Ca\ = (42.V“) GeV^. 


(51) 


Factorization implies that this contribution - up to the sign of Ca, and up to some 
constants, as explained at par. 2.2.2 - should be present also in the K decay. Here, 
however, also the current term is important [9], since, in this case, it receives contributions 
exclusively from the penguin amplitude, which is dominant over the tree amplitude. This 
fact is confirmed by the sizes of the proportionality coefficients of the current form factors 
for that decay [23]. Therefore, in order to calculate the current term, we insert into eqs. 
(|^ and (|^ the formulae of the electromagnetic nucleon form factors in the time-like 
region [9[ 123] , shown diagrammatically in fig. such formulae are given in the Appendix, 
while the proportionality coefficients of the form factors, Sq and Si, are listed in Table 


As regards the transition term, we assume the value (51) for Ca, with the two possible 
signs. The choice of the - sign yields 


Bf{B^ K^PP) = (2.96 ± 0.27) ■ 10"®, (52) 

which is considerably smaller than the experimental value (|^. The choice of the -|- sign 
for Ca yields even a lower value. Also data [ 2 ], [ 13 ] of the differential branching fraction 
are in disagreement with the factorization assumption, as shown by the dashed line in fig. 

therefore we conclude that this assumption is inadequate. The situation is somewhat 
analogous to the one illustrated in ref. [26] about the decays oi B to Kxo, to and 

to 7r®7r®, which are surprisingly enhanced with respect to the factorization assumption and 
demand both non-factorizing terms and re-scattering effects. 


3.2 Remarks on the Main Contributing Graphs 

After the result just shown above, a criticism of the assumptions made is in order. 
Therefore we analyze in detail the main graphs contributing to the decay considered. 



u 



Figure 1: Diagram of annihilation of a virtual 7 into pp 


which we illustrate in hg. 


Figs. 2a and 2b represent, respectively, the color-allowed and the penguin diagram. 
According to the former one, the meson is by no means correlated to the baryons, since 
the quark hadronization to the meson occurs just after the weak (very short-ranged) 
interaction; therefore factorization could be of course assumed for this diagram, which 
contributes to the transition term. But as told, in the K decay, this amplitude is sub¬ 
dominant with respect to the penguin one. This latter consists of several components 
contributing to the K and/or tt decay, described and represented in detail in refs. [27j 
and [7]: some of them may be regarded as factorizable, but not all. The evident failure of 
factorization leads us to conclude that non-factorizing diagrams play quite an important 
role in the K decay. In particular, we point out a contribution fed by the quark subprocess 
b —)■ CCS. It may give rise to a pp pair through re-scattering from pair formation, 

which can occur either by cc fragmentation, or via intermediate states like 
2 cAc[Z]; this causes an anomalously large enhancement of the decay rate with respect to 
the naive factorization [7]. In the vr decay, the corresponding amplitude is suppressed by 
a factor |V/(i/V/s| with respect to the K decay. Since the current term interferes only 
weakly with the transition term[S], this process gives a negligibly small contribution to 
the differential decay width, as anticipated in subsect. 3.1. On the contrary, it might 
explain the prevalent role of the current term in the K decay, as we shall establish in 
subsect. 3.3. 

The failure of factorization is also conhrmed by the numerical results presented by 
Chua et a/. [S], since the current contribution depends crucially on the effective number 
of colors, Nc (see also ref. CZD- On the contrary, the £t to the tt decay, dominated by 
the tree diagram, is stable versus Np., this is reflected also on the transition term of the 
K decay [9], for which those authors assume the same parametrization as for tt. The £t 
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Figure 2: Main diagrams in the decays —)■ K^pp 


to —)■ 71^Ap, which depends essentially on a factorizable penguin diagram |27j. is also 

reasonably stable versus Nc. 

Aside from that, as told, the long-range FSI might be not so negligible, somewhat 
analogously to the case of B decays to mesons[26], although on a smaller scale. This 
may contribute to the failure of factorization and may give rise to further effects, to be 
discussed in sect. 4. 

3.3 An Alternative Parametrization 

All of the above considerations suggest to relax the factorization assumption. Indeed, we 
have to take into account both non-factorizable terms and FSI; these latter are typically 
non-perturbative and therefore cannot be deduced from QCD principles. To this end, we 
regard Ca, Cy;,, Cp, Sq and Si as free parameters, which we determine by minimizing 
the in a fit to data, that is, as told, the differential branching fraction [21 [15] and the 
asymmetries App and Aqp^- In order to avoid strong correlations between the various 
parameters, we set Cvb = Cp = 0, like Chua et a/. [D]. The parameters of the fit are shown 
in Table [T| where they are compared to those of factorization. Furthermore the fit to the 
rupp distribution is exhibited in fig. [^ continuous line. 
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Figure 3: K^pp decay: differential branching fraction vs rupp, charge averaged 

case. LHCb and Belle data are represented, respectively, by rhombuses and triangles. 
The continuous (dashed) line refers to the best £t (factorization assumption). 


Table 1: Parameters of the theoretical function. FI refers to factorization, F2 shows the 
results of the best £t. Ca and Cys are expressed in GeV^, Cp m. Geld®, S'o and Si in 
Gel/. 


Fit 

Ca 

Cv5 

Cp 

So 


FI 

-42.1 

- 

- 

0.36 

0.36 

F2 


- 

- 

1 Q7 q+ 48.2 
lO 1 .0_]^26.9 

1 or+O.OS 


3.4 Results and Model Predictions 

Table shows the experimental values of the branching fraction and of the CP and 
FB asymmetries and the theoretical results about such quantities, both according to 
the factorization assumption and to our best £t. Furthermore, hgs. and represent. 
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respectively, the differential branching fractions of the and K~ decays separately and 
the differential CP-asymmetry, Acp, eq. (49), which could be measured in the future. 

Two comments are in order about the results given in Table and in Table First of 
all, the parameters Ca and S'o, as determined by our best fit (see line F2 of Table [^, are 
affected by large, asymmetric errors; however, while the value of Ca is compatible with 
factorization (compare with line FI), So is much greater, which suggests that the current 
term derives a large contribution from the effects described in subsect. 3.2. As regards 
Table 1^ the theoretical value of Acp is affected by a large error, while the one of Aps 
is quite small. This is connected to the experimental errors and to the definition itself of 
these asymmetries. Acp is proportional to the difference between two large and almost 
equal quantities and is known experimentally with a relative error of 141%; therefore, its 
dependence on the fit parameters is quite mild. The opposite happens for App, whose 
experimental value is not so close to the theoretical one and is affected by a quite small 
relative error [2]. 


x10‘® 



Figure 4: Differential branching fractions vs rripp for —)■ K^pp decays separately: the 

continuous (dashed) line refers to the decay. 
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Table 2: Branching fraction and CP and FB asymmetries: experimental data (Ex) are 
compared to results of the factorization assumption (FI) and of our £t (F2). Bf is 
multiplied by 10®, Acp and Afb 10^. Here the experimental errors of the asymmetries 
are summed quadratically (see eqs. (j^ and (j^). 


Fit 

Bf 

Acp 

Afb 

Ex 

5.9 ±0.5 

-2.2 ±3.2 

37.0 ±2.4 

FI 

2.96 ±0.27 

o n +0-24 

O-U -0.29 

7.6 ±0.1 

F2 

6 21 

-0.49 

11.7 

11.41?:' 


4 Discussion and Conclusions 


Here we comment on our main results. 

a) As regards the current contribution, we have adopted the parametrization suggested 
by the pole model|23], like ref. [9|. As already observed in subsect. 3.4, by comparing 
the lines LI and L2 of Table we conclude that, according to our £t, the current term 
is much greater than expected from factorization. Moreover, we see from line F2 that 
|S'o| >> I S'! I, therefore the term corresponding to the angular momentum J = 0 prevails 
neatly over the one with J = 1, in agreement with the low-mass character of the pp peak. 
However, the term with J = 1 gives an essential contribution to the £t to the Dalitz plot 
distribution, as expected from C and P violation in weak decays 

b) The FB asymmetry {Afb) - confirmed by the Dalitz plot distribution [21 HZ] - is 
caused by the linear terms in cos dp, which appear in the differential decay width 
precisely. 


- in the Jacobian (36), owing to the Lorentz boost from the pp frame to the B rest 
frame; 

- in the interference term between the pseudoscalar and the vector amplitude, whose 


coefficient is Dis, eq. (31). 

In our model - for which, incidentally, the second cause of asymmetry shown above is 
ruled out, owing to the absence of the pseudoscalar amplitude - Afb is positive for the K 
decay, not only as a result of our £t, but also according to the factorization assumption, 
as shown in TableIn this connection, it is worth observing that, in the model proposed, 
the sign of Afb is not a priori determined; it depends on the moduli and phases of the 


amplitudes involved in the decay, according to the (3 coefficients (20) to (23). Furthermore, 


this is not typical of the K decay. Indeed, also in the case of tt decay, the parametrization 
by ref. |9] yields a positive Afb-, whereas experiments hnd a negative value [TTl IT^ [2]. 
Therefore, our picture seems to contradict the pole model conclusions and naive quark 
model expectations [7|. Rather, as observed in subsect. 3.1, the negative FB asymmetry 
of the TT decay could be explained by means of the color-suppressed tree diagram. 

c) In the K decay, the gap between the experimental value oi Afb and the result of our 
fit (see Table 1^ could be filled by the first term of eq. (22), concerning the electric form 
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Figure 5: Differential CP asymmetry vs rupp 


factor of the proton, neglected also in ref. |9], which is linear in cos 6p. This explanation 
of the large asymmetry could be alternative to the one proposed by other ant hors [TT], 
d) Lastly, the CP asymmetry {Acp) is found to be positive according to our 
parametrization, which does not differ so much from the prediction of the factorization 
assumption (see also ref. m) ; on the contrary, experiments seem to indicate a negative 
value, although compatible with a positive asymmetry within errors. Therefore, if a more 
accurate determination of this observable could be desirable, also a deeper analysis of the 
model is needed. Indeed, although less important than in meson B decays, the FSI could 
produce a relative phase between the current term and the transition amplitude, so as to 
change the sign of Acp |26l |28]. In other words, the difference of sign between data and 
model predictions does not necessarily imply new physics. 
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Appendix 

Here we give the numerical values of the CKM matrix elements and of the Wilson 
coefficients involved in the decays considered, as well as the expressions of the matrix 
elements that appear in the expressions of the current and of the transition terms, 
according to refs. [91123]. 

CKM Matrix Elements 

We give the CKM matrix elements [23] which enter the expressions of the decay amplitudes 
for B~ — K~{n~)pp: 


14i, = AX^p-iri), V:, = \ (A. 1) 

V- = -AA^ lC, = l-iA^ V;, = AX^{l-p + tp), (A. 2) 

with A = 0.814, A = 0.226, p = 0.117 and rj = 0.353. For the CP-conjugated decays one 
has to take the complex conjugated elements. 

Wilson Coefficients 


For the —)■ K^pp decay, we have]^] 

ai = 1.05, 02 = 0.02, 

04 = -(3.873 + 1.21^)■10-^ 
06 = -(5.553 + 1.2H)-10"^ 


03 = (72.7- 0.3i) ■ 10■^ 
ag = -(66 + 0.3*) ■ 10-^ 
09 = -(92.6 + 2.7*) ■ 10"^ 


As for the B^ —)■ n^pp decay, the coefficients read as 


oi = 1.05, 02 = 0.02, 

04 = -(3.757 + 1.083*) ■ 10-^ 
06 = -(5.447 + 1.083*) ■ 10“^ 


03 = (73 + 0.3*) ■ 10■^ 

05 = (-66 + 0.3*) ■ 10-^ 
09 = -(92.4 + 2.5*) ■ 10-h 


(A. 3) 


(A. 4) 


Non-perturbative Matrix Elements 

A) The matrix element {K~\L^^\B~), which appears in the current term ([^, reads as[9] 

{K-\L^jB-) = k^^rFo + BF^, (A. 5) 
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with 


and 


Fn = 


k = Pp + pp, I = 2pB - (1 + r)k, 
So ^ S, 


m 


r = 


B 


m 


K 




(A. 6 ) 


1 — UiT + U2 t‘^ ’ 


Fi = 


(1 — r)(l — wr) ’ 


^ (A- 


m 


V 


In order to get the expression of (tt \L^i\B ), we have to replace by in r. Moreover, 


the numerical values of the parameters which appear in eqs. (A. 7) are 

mv = 5.42 (5.32) GeV, Ui = 0.70 (0.76), U 2 = 0.27 (0.28), w = 0.43 (0.48), (A. 8 ) 

the values outside (inside) parentheses referring to the Kij) decay. Lastly, Sq and Ai, 
which have the dimension of an energy, are treated as free parameters and determined 
according to the two possible choices described in the text: see Table 

B) Now we give the expressions of other matrix elements involved in the current terms 
in K decay [9]: 


{K\L,,\B) = 


- rr 


mb - ms 


0 , 


{pp|L;,| 0) = u(pJ|<l>-7„ + i>‘ ■ ■- 


2 m, 


- i.9Al^llb + 75X- k^)]v{pp) 


2mp 


{pp\Rf,\0) = u{pp)[^^jp + + 75; 


{pp\kss\Q) = {—g\ + 


m.s 


Ampm, 


2mr 


-h\)u{Pp)l5v{pp). 


2mp 


kp)]v{pp) 


(A. 9) 
(A. 10) 
(A. 11) 
(A. 12 ) 


Here the first and last formulae are deduced from equations of motion; in particular, 
the hrst one exploits also eq. (A. 5), contracted with k^, and the relation k ■ I = 0 . 
Moreover 


^ Pp Pp) 


(A. 13) 


$1 

$3 


mP 


m 

Xi 

X3 

X5 

yi 


= + <l>2 = 3(5y-G^), 


c mP 

Vi)E 


2 G;„ 4' = G^, 


n=\ 


Xn 

^n+ 1 ’ 


Gl, = -Mt) Y. 


n=l 


Un 

^n+ 1 ’ 


Sv = -A(i)|, 


[ln{t/Al)]\ Ao = 0.3 GeV, 7 = -2.148, 
420.96 GeV\ X 2 = -10485.50 GeH®, 
106390.97 Gel/®, X 4 = -433916.61 Gel/^°, 
613780.15 Gel/^2 


(A. 14) 
(A. 15) 
(A. 16) 

(A. 17) 
(A. 18) 


= 236.69 GeV^ 


= -579.51 Gel/®, = -52.42 Gel/^ (A. 19) 
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and 


9\ — f A ~\~ dA ~\~ g\ — 2dA + 3sa, 

Qa = ~fA + dA + sa, 9\ = fA + l/Sc^A, (A. 20) 

AA-mA (A. 21) 

n=l n=l 

wi = 399 Ge1/^ W2 = -1055 Gel/®, (A. 22) 

= 65.93 GeV^, Z 2 = -1055 Gel/®, Z 3 = 333.06 Gel^^ (A. 23) 

Lastly, we have set 

^ — AiTfl^ 

= 0 and h\ = - —9a, i = 1 to 4. (A. 24) 

t 


The corresponding expressions for the tt decay are obtained by replacing mx by m^r, 
by irid, by Sy and g\ by s^. 

C) The matrix elements involved in the transition term @ areP] 


(/L|L„s|0){pp|S'“'’|0) = -fKmb[FAu{pp)^K'y5v{Pp +Fpu{Pph5v{pp)], (A. 25) 
(/L|L„s|0)(pp|S'“'’|0) = -fKmb[Fv5u{Pp)^Kv{pp +Fpu{pp)v{pp)]. (A. 26) 


Here fx = 0.158 GeH and 


Fa = 


Ca 
F ’ 


Fv5 — 


Cv5 
F ’ 



(A. 27) 


Ca, Cv 5 and Cp being assnmed, again, as free parameters and determined as described 
in the text. 

In the case of the tt decay, one has to replace fx by fT, = 0.133 GeV. 
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